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Abstract 

^ We consider a ballistic Mach-Zehnder interferometer for electrons propagating chirally in 

one dimension (such as in an integer Quantum Hall effect edge channel). In such a system, 
dephasing occurs when the finite range of the interaction potential is taken into account. Using 
the tools of bosonization, we discuss the decay of coherence as a function of propagation distance 
and energy. We supplement the exact solution by a semiclassical approach that is physically 

S transparent and is exact at high energies. In particular, we study in more detail the recently 
predicted universal power-law decay of the coherence at high energies, where the exponent does 
>-^j not depend on the interaction strength. In addition, we compare against Keldysh perturbation 

theory, which works well for small interaction strength at short propagation distances. 
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1 Introduction 

The loss of quantum mechanical phase coherence by a fluctuating environment plays an essential 
role in many different branches of modern physics [T] . It governs the transition from the quantum to 
the classical world [SJ [5] and occurs as an unavoidable consequence of any measurement process [1] . 
It introduces the dephasing time as the time-scale during which it is possible to observe quantum 
coherent dynamics in qubits and other quantum objects. Furthermore, the interference contrast 
in man-made interferometers is determined by decoherence. This includes setups for electrons in 
semiconductors or metals, for single photons, neutrons, neutral atoms and larger objects. In each 
of those examples various different sources of fluctuations contribute, e.g. thermal radiation, or the 
interaction with phonons and other particles. 

The paradigmatic setup that is treated in much of the literature on decoherence is the following: 
A "small" quantum system with few degrees of freedom or even a finite-dimensional Hilbert space 
(atom, spin, qubit, single particle) interacts with an equilibrium environment that comprises an 
infinite number of degrees of freedom, i.e. a "bath". In the majority of cases, the model of the 
environment is restricted even further, to consist of a collection of harmonic oscillators (e.g. photon 
or phonon modes). Quantum dissipative systems of that kind already offer a rich phenomenology, 
including exponential or power-law decay of the coherence in time, as well as dissipative phase 
transitions at strong coupling [5] [T] . 

However, there are several situations in which one is forced to go beyond that class of models. 
We mention three of the most important examples, where the nature of the environment has to 
be reconsidered, (i) The environment may be driven out of equilibrium (as in the interaction with 
electrical currents or a photon stream emanating from a laser) , which is particularly important for 
measurement setups. As a consequence, extra fluctuations are introduced and there is no simple 
relation any more between the dissipative response and the fluctuations, (ii) The environment 
may be different from a bath of harmonic oscillators, such that the resulting fluctuations are not 
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Gaussian. This comprises examples like spin baths, fluctuators in solids, or potential fluctuations 
being produced by discrete particles, such as electrons. The peculiar features of non-Gaussian 
environments may be observed at strong coupling, and they include oscillations of the interference 
contrast with time [SlU]. (hi) It may be hard to cleanly separate (conceptually and technically) the 
environment from the system. This is the case for interacting many-particle systems, where one is 
interested in the effects on a single particle that is part of the whole, being indistinguishable from 
the other particles. 

In this paper, we will address the third case, which is important in a variety of electronic systems, 
ranging from dephasing in disordered conductors to electronic interferometers. Here we will focus 
on interacting electrons moving inside a ballistic interferometer, although we also would like to 
mention our recent related work on the importance of Pauli blocking for dephasing in disordered 
electron systems [5J [TU] • More specifically, we will discuss the loss of coherence in an interacting 
one-dimensional chiral fermion system. Besides representing an exactly solvable model system, 
where one is not restricted to perturbation theory in the interaction, this also has become relevant 
for experiments recently. In a series of experimental studies [TTJ Q21 H3J HH [H] , first initiated at the 
Weizmann institute, an electronic Mach-Zehnder interferometer has been realized, employing edge 
channels in the integer quantum Hall effect. The interference contrast as a function of voltage and 
temperature has been analyzed, and only a fraction of the features have been explained by now. 

On the theoretical side, dephasing in such a setup has been discussed both for dephasing by 
external fluctuations [TU El HHJ EH 123 (2D EH IS] (such as phonons, defect fluctuators or Nyquist 
noise from external gates, or "dephasing terminal" reservoirs), as well as by the intrinsic electron- 
electron interaction |2"2"1 1231 R 1231 1251 |251 |2T] . 

Electron-electron interactions in one-dimensional systems are usually described within the Lut- 
tingcr liquid framework. Dephasing of electrons in Luttingcr liquids is interesting as an example of 
a non-perturbative, strongly correlated model system, and has been studied already in a number 
of works 28, 29 3U1 [31] |321 [33J • In contrast, the situation for (spinless) chiral interacting fermion 
systems, such as edge states in the integer quantum Hall effect (QHE), seems to be straightfoward. 
Within the commonly discussed "g-ology" framework, one assumes pointlike interactions. In that 
case, an interacting chiral model is nothing but a Fermi gas with a renormalized velocity. No 
further interaction effects are present, in particular there is no dephasing. However, for interferom- 
eters in nonequilibrium (at finite bias voltage) or at finite temperatures, one may still encounter 
interesting physics, provided the finite range of the interaction potentialis accounted for. This was 
discussed only recently [23], although there are some nice and thorough earlier studies [3U [35] of 
the momentum-resolved density of states in such systems. 

In the present paper, we build on the recent analysis of Chalker, Gefen, and Veillette (25]. These 
authors modeled the interacting Mach-Zehnder interferometer as consisting of two one-dimensional 
chiral interacting fermion systems, corresponding to the two channels or arms of the interferome- 
ter. The two channels are tunnel-coupled weakly at two locations, representing the quantum point 
contacts (i.e. the beam splitters) of the experimental setup. By staying within the regime of weak 
tunnel coupling (low transmission), one is able to express the current in lowest order perturbation 
theory with regard to the tunneling. The result is a formula for the current (and, subsequently, 
the interference contrast or "visibility" ) that only involves the Green's functions of the interacting 
channels in the absence of tunneling. These Green's functions can be obtained using the tools 
of bosonization. In order to obtain nontrivial results, it is necessary to go beyond the commonly 
employed assumption of pointlikc interactions. Only when treating the full dependence of the in- 
teraction potential on the distance of the electrons, interaction effects beyond a simple velocity 
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renormalization are observed. The main results of their study are that at low voltages and tem- 
peratures the interference contrast becomes perfect, while the suppression of contrast at increasing 
voltages and temperatures depends on the details of the interaction potential. 

In our work, we will first review the general expression for the current that consists of two parts, 
a flux-independent term and the interference contribution that will be suppressed by interaction- 
induced decoherence (section |2| . In contrast to [25] , we formulate the answer in terms of the 
Green's function in energy-position space. This has the advantage of corresponding directly to 
the contribution of an electron at energy e that travels a distance x inside one of the arms of the 
interferometer. Next, we review the model Hamiltonian and the solution by bosonization (section 
|3J). We discuss the general features of the Green's function \G > (e,x)\ that has been obtained by 
numerical evaluation of the exact bosonization expressions, and study the influence of the coupling 
constant. At low energies, the decay with propagation distance x is weak. It becomes faster when 
the energy rises above the characteristic energy that characterizes the finite range of the interaction. 
Finally, at large energies, the decay becomes independent of energy. 

The latter asymptotic regime is then the subject of a more detailed analysis using a physically 
transparent semiclassical picture (section [4]) . This picture of a single electron interacting with the 
potential fluctuations produced by the other electrons has been exploited by us in a recent short 
paper |36j to discuss dephasing in an interacting chiral fermion system. It had been introduced 
earlier to deal with dephasing of ballistically propagating electrons in contact with an arbitrary 
quantum environment |17l 121] , and has also been suggested independently in the context of two 
interacting Luttinger liquids [3U]. Here, we provide more details of the calculation and an extended 
discussion of the fluctuation spectrum that is seen by the moving electron in its frame of reference. 
One of the main results is that at high energies there is a "universal" power-law decay \G > (e, x) \ oc 
1/j of the electron's coherence, with an exponent independent of interaction strength. We also 
analyze the situation at finite temperature, where one has to discuss the transformation of the 
lab- frame temperature into an effective temperature in the co-moving frame of reference. 

Although in principle, for this particular problem, the bosonization solution is fully sufficient, we 
conclude our analysis with a section on perturbation theory (section [5j. This is done in anticipation 
of going to different setups where an exact solution is no longer possible. Even though we will only 
discuss the equilibrium Green's function of the interacting system, we employ the nonequilibrium 
(Keldysh) diagrammatic technique, to provide for a straightforward extension to situations where 
this is needed. We calculate the self-energy up to second order in the interaction. This includes a 
diagram that describes the decay by emission of a plasmonic excitation, which is however partially 
cancelled at low energies by an exchange-type diagram. We will show that at short propagation 
distances and for modest coupling strengths, the Keldysh result provides a good approximation to 
the exact solution for |G > (e,a;)|, even though its structure in (e, fc)-space is qualitatively different 
from the bosonization result. 



2 The electronic Mach-Zehnder Interferometer 

The electronic Mach-Zehnder interferometer is one of the simplest model systems where the inter- 
play of quantum mechanical coherence and many-body effects can be studied, both in theory and 
experiment. 

To set up the description, we first neglect interactions and imagine a single electron moving 
through the interferometer. The interferometer itself is described as two parallel one-dimensional 
channels in which electrons propagate into the same direction (see Figflk) . At two tunnel contacts 
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(i.e. quantum point contacts (QPC's) in the experimental realization), these channels are coupled 
by tunneling amplitudes tA and tg. Further below, we will assume these tunneling probabilities to 
be small perturbations, coupling lead 1 (left channel) and lead 2 (right channel). Furthermore, a 
magnetic flux is enclosed by the interferometer, which leads to an additional Aharonov-Bohm phase 

4>. 

In the experiment the current / through the interferometer, i.e. the current between the two 
leads, measured at the output port, is the quantity of interest (see Figjl]). It contains two types of 
contributions: one flux- independent constant term and one interference term that depends oncos(0). 
The contrast of the interference fringes observed in I((f>) = Iq + I C oh(<f>) can be quantified via the 
so-called visibility 

^max -*min /-, \ 

vi = j TJ — (1) 



L mill 



(^min) are the maximum (minimum) current as a function of flux. This definition is 
chosen that the visibility is equal to one for perfect interference contrast. This can be used as a 
direct measure for the coherence of the system. The coherence can be destroyed by the influence 
of an external bath as well as by internal interactions like the Coulomb interaction between the 
electrons inside the interferometer. 

Treating the interferometer as a many-body system yields expressions for the current through 
the interferometer, which are not as obvious in a physical sense as in the single particle picture. 
Therefore the goal of this section is to formulate the quantities of interest in the physically most 
intuitive and transparent way. 

2.1 Current 

The observable of interest in the present setup is the current through the interferometer (see 
Fig{TJi)due to a finite bias voltage between the two leads, i.e. ^i — ^2 ^ 0. Dealing with the 
electron-electron interaction exactly using the bosonization technique has the disadvantage that we 
have to treat the tunneling between the channels in perturbation theory. In the following we define 
the quantum point contacts A and B to be at the positions x^ = and x^ — Xj, respectively 
(where j — 1 , 2 is the channel index) . Then the tunneling Hamiltonian is given by 

H T := i A fct(0)¥ a (Q) + t B %l(xi)% 2 (x 2 ) + h.c. (2) 
The current into channel 1 is (q e < 0): 

I=^j t Nx. (3) 
Thus the Hciscnbcrg equation of motion yields 

7=- ge i[^ 1 ,ff + ff T ] = -q e i[^# t 1 (0)# 2 (0)+i B * t 1 (x 1 )# 2 (x 2 )] +h.c . (4) 

Now we change to the interaction picture with respect to Hq (where Hq denotes the interacting 
electron Hamiltonian in the absence of tunneling, which we define in section p|, setting AH (t) = 
e iH t A e -iH t _ -yy e are interested in the steady-state current through the interferometer, which we 
obtain as a Kubo-type expression, in linear response with respect to the tunneling Hamiltonian, at 
arbitrary bias voltages: 

I=)f di([lff o (0), -ffffo (*)])• (5) 
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Figure 1: (a) Scheme of the interferometer setup. The two channels 1 and 2 of length and 
the corresponding chemical potentials are indicated. The electrons can tunnel at QPCs A and 
B, with tunnel amplitudes t& and ts- By tuning the magnetic flux $ through the interferometer, 
one observes an interference pattern /(</>). The solid orange lines denote the current through the 
interferometer. 

(b) A single electron propagating at high energies feels a fluctuating quantum potential V(t), due 
to the interaction with the density fluctuations in the sea of other electrons. We show a density 
plot of the potential, which is produced by the electronic density fluctuations (plasmons) in the 
channel. The plasmons are moving with a renormalized velocity v (see main text) while the high- 
energy electron moves with the bare Fermi velocity vp. It picks up a random phase, which leads to 
dephasing. This is the picture underlying the semiclassical approach in Section [4] 
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2.2 Evaluation of the current 

Starting from Eq. ([5| and plugging in the apropriate definitions an explicit expression for the 
current can be found. It will be convenient to express the current in terms of the unperturbed 
Green's functions of the chiral electron liquids in the two channels. These Green's functions will 
therefore be the primary object of our subsequent discussions. 

2.2.1 Green's functions 

In particular, we will analyze the particle- and hole- propagators 

G>(x,t) = -i(i>(x,t)¥(Q,0)), (6) 

G<(x,t) = +i^¥(0,0)^(x,t)\ (7) 

where we omit the channel index for brevity, unless needed for clarity. The Fourier transforms are 
obtained as G > (k,uj) = J dx J dte~ l ( kx ~ wt 'G > (x,t). In addition, we will need the spectral density 

A(k,u) = --Im[G R (k,ui)], (8) 

where G R (x,t) = -iO(t) (^$(x,t),& (0,0)}^. The energy-dependent tunneling density of states 
follows as 

/oo 
dkA((J,k). (9) 
-oo 

2.2.2 Flux-independent part of the current 

The flux-independent part of the current is found using Eq. @, @ and 

Io = qe(\t A \ 2 + \t B \ 2 ) f°° dt [G>(0,-t)G<(0,t)-G<(0,-t)G>(0,t)], (10) 

J — OO 

which we rewrite by going to the frequency domain 

J = q e (\t A \ 2 + \t B \ 2 ) J (dw) [G>(Q,oj)G<(Q,u;) - G<(Q } u;)G>(0, w)] , (11) 
where J (dui) . . . = J du)/2ir .... Using the identities 



Gf(x = 0,u) = J {dk)Gf(k,w) = -i[l-f j (u)]-u j (u) (12) 

Gf(x = 0,w) - i/j-H-i/^w), (13) 

where fj{oj) denotes the Fermi function. We can reformulate the expression as a function of the 
tunneling density of states ^(w), 



Jo = ?e {\t A \ 2 + |<b| 2 ) / {duS)vx{u)v 2 {u) 



/2(w)[l-/iM]-/i(w)[l-/a(«)l 
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which finally yields the most intuitive form describing the sum of tunneling currents at two point- like 
locations: 

h = Qe (\t A \ 2 + \t B \ 2 ) J (dw)i^(w)i^(w) [/ 2 H - AM] . (14) 

In particular, these expressions show that the flux-independent part of the current only depends 
on the tunneling density of states. It is independent of the length of the interferometer arms. This 
is to be expected, as that part of the current is insensitive to the electrons' coherence, and therefore 
the decay of coherence as a function of propagation distance will not enter here. 

2.2.3 Coherent part 

The Mach-Zehnder setup is intended to investigate the coherence of the electron system and there- 
fore the most interesting quantity is the coherent part of the current, which we define to be the 
flux-dependent contribution. Again, using Eq. @ and Q it yields : 

/coh(<£) = Qc J (duj) [(t A t* B )e- 1 * ■ G>(lu, Xx)G<{u, -x 2 ) - (t* A t B )e*+ ■ G<(lu, -x 1 )G>(cj,x 2 )] + c.c. . 

(15) 

At T = 0, in a situation where the particle current flows from channel 2 to 1, only the first term 
(and its complex conjugate) contributes. 

It might be helpful to see how the structure of this term G^G 2 can be understood in an intuitive, 
if slightly imprecise, way, that also relates to our subsequent semiclassical discussion. When the full 
beam in channel 2 impinges onto the first beam-splitter A, we obtain a superposition between two 
many-particle states: With an amplitude near unity, nothing happens (no tunneling takes place), 
and we denote this state as |0). There is a small chance (of amplitude t A ) for a particle to tunnel 
through A into channel 1, leaving behind a hole in channel 2. As time passes, the second part 
acquires an amplitude (relative to the first) that is given by the product of propagation amplitudes 
for the electron (ipi) and the hole (ip 2 )i resulting in: 

|0) + UV>iV4 |lp,2fc) . (16) 

Including the action of the second beam splitter B, and the Aharonov-Bohm phase, the total 
probability to detect an extra electron in the output port (channel 1) is therefore 

\t B e** + t A ih.il>2\ 2 , (17) 
which gives rise to the interference term 

t A t* B e-^^* 2 +c.c. (18) 

Averaging the amplitudes over phase fluctuations induced by the interaction, we arrive at the 
propagators, replacing (i/'i) by the particle propagator Gf, and (tp^) by the hole propagator G 2 ■ 
The full analysis keeps track of energy conservation 

Thus, in the many-body picture, the observation of an interference term in the current is seen 
to depend both on the passage of an electron through channel 2 (G 2 ) as well as the coherent 
propagation of the corresponding hole, of the same energy lu, in channel 1 (Gf). This issue has been 
discussed before, both for the Mach-Zehnder interferometer and weak localization [TT1 |2"T1 ITTTl 15] . 
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2.3 Visibility 

In the Mach-Zehnder setup, the so called visibility is used as a measure of the coherence of the 
system. There are different definitions (experimentally, often the differential visibility is employed). 
However we will define the visibility in terms of the total current, as 

. . _ max0/ co/t (</)) 

^max -^min /1 r\\ 

v I {V,l)= = - — — . (19) 



max 



The bias voltage is defined as [i\ — /i 2 = q e V and we set V > 0. We will also focus on zero 
temperature, as this seems to be the most interesting case. Now the visibility can be written in a 
compact form (here shown for T = 0): 



2\t A t 

vi - 



J^ Vl du G>(u,x x ) -G<(lo~ \g e \ V, -x 2 ) 



M 2 + \t B \ 2 du; u(u) ■ v{u - \q e \ V) 



(V > 0). (20) 



Note that the channel indices of the Green's functions (GFs) are omitted, as in this formula the GFs 
are defined with respect to a fixed density (and potential Uj = 0) and all the explicit dependence 
on the bias voltage is shifted to the GF arguments. 

Motivated by the structure of Eq.[l5]in the following we will focus the attention on analyzing the 
function ' (e, x) in three different ways. First we will apply the bosonization technique, i.e. we will 
include the intrinsic interaction in all orders. The main disadvantage of the bosonized picture is that 
we are looking at single particles tunneling between channels, while we are phrasing the description 
in terms of collective, bosonic excitations which prohibit a simple physical picture of the process of 
dephasing. Therefore, as an alternative point of view, we will discuss a transparent semiclassical 
model for electrons propagating high above the Fermi energy, subject to the fluctuations produced 
by the rest of the electrons. Finally, we will complement this analysis by studying the system in 
Keldysh perturbation theory for the fermions, up to second order in the coupling strength. 



3 Solution by bosonization 

In this section we recall how to obtain the Green's functions needed here via bosonization [571 155] , 
Experienced readers may skip the section, and can refer to it later regarding the notation. 



3.1 Hamiltonian and formal solution 



We start from interacting right-moving chiral fermions in the two channels j = 1,2, employing a 
linearised dispersion relation (we set H — 1): 



.7=1,2 



J2( U 1 + V Fk)c] k C^ k + tfintj 
.fc>0 



(21) 



where vf denotes the Fermi velocity and Uj is a constant which fixes the chemical potential of the 
channel. The particle operators for the chiral electrons are 



Cj.k ■ 



(22) 



fc>0 
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where L describes the size of the artificial normalization volume (with L — > oo in the end) . We also 
introduce the density fluctuations within each channel: 

p j (x) = ^ j (x)^ j (x)-j5 j (23) 

such that the Fourier components are given by 

= tE^^ ( 24 ) 



Pq 



fc>0 



with f>-qj = fiq j. The average density pj enters pj = Uj + 2irvFp~j. 

As we take the two interferometer channels to be spatially separated, we only have to take care 
of intrachanncl interactions. Transforming the interaction Hamiltonian 

A^-I/*/**}(.)*jMi7<.-^*,M*,M w 

into momentum space yields 

E u A+ q ,A'- q ,j tk ''i tk >i' ( 2? ) 



fc>0,fc'>0,<j 

where U q — J dxe~ lqx U{x) are the Fourier components of the interaction potential. 

We construct bosonic operators from the Fourier components of the density in the standard 
way. As we are only dealing with chiral electrons, we only need to consider q > in the following: 

/ 27r \ 1 /2 /27r\ 1/2 

(<? > o) Ki = [jq) p 9 j = [j^) P-™ ■ ( 28 ) 

These operators fulfill 

ikvfyJ = ^"W • ( 29 ) 

As usual, the main advantage of bosonization consists in being able to write the kinetic part of the 
Hamiltonian as a quadratic form in boson operators: 



.3=1,2 



<?>0 



+ const . (30) 



The interaction part of the Hamiltonian reads: 



L 

q>0 
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Thus, the Hamiltonian is already in diagonal form, 



i=i,2 



.g>0 



with the plasmonic dispersion relation of an interacting chiral ID electron system: 



uj(q) = v F q 



2nvp 



(32) 



(33) 



For the following discussions, we introduce the dimensionless coupling constant a = ' 7 2 ^^°- ) , where 
a g] — 1, oo [. The renormalized plasmon velocity at small wavenumbers is v = vp(l + a). Negative 
values of the coupling constant are related to attractive interations, positive values to repulsion (at 
small wavenumbers). For a — > —1 the plasmon velocity tends to zero, v — > 0. For a < — 1 the 
system is unstable, i.e. formally uu(q) < for q > 0. 

The final step is to express the single-particle operators using bosonic fields [38l [37] (assuming 
L — > oo): 



Jk F x -i<ij(x) 



q>0 



(34) 



The ultraviolet cutoff length a is sent to zero at the end of the calculation. The Klein operator Fj 
annihilates a fcrmion in a spatially homogeneous way, with the following commutation relations: 



{Fu FJ} = 6 id ; {Fi, Fj] = {fIfH = 



(35) 



The explicit time dependence of the Klein operators is obtained by using the Heisenberg equation of 
motion, i.e. dFj/dt = —i[Fj,H] — —i/ijFj Fj(t) — e~ lfJ - jt Fj(0). In the end taking into account 
the relation kp j — we arrive at 



* j (x) = F j 



'27ra 



(36) 



3.2 Green's function from bosonization 



Now we are able to evaluate the Green's function defined above explicitly using the bosonized single 
particle operators The calculation is done quickly using the fact, that the Hamiltonian in terms 
of the bosonic operators is quadratic, i.e. the field Q[b, tf] can be treated like a Gaussian (quantum) 
variable, resulting in 



G>(x,t) 



-i/j,(t— x/vf) 



2wa 



exp[($(M)$(0,0))- (*(0,0) 



(37) 



and analogously for G < . By factoring off the non-interacting Green's function Gq^, we can write 
G >/< {x, t) = e-W*-*/"*] • G> f< (x, t) ■ exp[S* fl (a;, t) T iSi(x, t)] , (38) 

1 1 



where 



G> /< (x > t) = 



2f3v F sinh[^- (a; - v F t ±i0+)] 



(39) 
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with /3 = A and fcs = 1. All the effects of the interaction now are included in the exponent where 
we have to subtract the non-interacting contribution: 



S R = 



— {coth[-^] [cos[w q t - qx] - 1] -S R (u q -> qv F )} 

o i j , 



/*°° dq 

Si = I — {s'm[uj q t - qx}-Si(uj q qv F )}. 
Jo 1 " -v ' 



(40) 



(41) 



3.3 Discussion: Green's function in space and time 

Here we discuss the Green's function as a function of space and time. A more detailed discussion 
can be found in [55]. The absolute value of the Green's function G > (x, t) is shown in Fig[2] at zero 
temperature (to which we restrict our discussion). In the following all the numerical evaluations are 
performed using a generic smooth interaction potential U q = Uoe~" 9 ^ qc ^ . We note that all the 
qualitative results are valid for potentials which are finite at zero momentum (Uq 0) and which 
are cut off beyond some momentum scale. Those assumptions are not restrictive and for example 
are fulfilled for a Coulomb potential with screening in a quasi one-dimensional channel of finite 
width. 

In Fig(2j we observe as the main feature that the Green's function splits into two parts during 
its propagation. One of those propagates with the bare Fermi velocity v F and represents the 
unperturbed Green's function, i.e. the high energy part. For increasing time its weight decreases, 
i.e. the amplitude of the bare electron to arrive at x without being scattered decreases. The other 
peak represents the low energy part, stemming from energies below e — fi ~ vpq c . It moves with 
the renormalized velocity v. 

We can obtain this structure of G > (x,t) from a crude approximation. Namely for T = in a 
first approximation we can cut the momentum integral due to the fact that for g> q c the integrand 
vanishes, i.e. 



S[x, t] 



iSi 



dq 



h <1 

The integrals are known and yield 



e~ q/qc [exp[-iq(vt - x)] - 1] 



dq 



e~ q/q " [exp[-iq(v F t - x)] - 1] 



S[x, t] « In 



x — vpt + iq c 
x — vt + iq^ 1 



Therefore the structure of the Green's function is given by 



(42) 



1 



x — vpt + i0+ 



x — vpt + iq c 
x — vt + iqc 1 



(43) 



displaying both the 6 peak at x — v F t and the broadened peak at x = vt, whose width is set by 
q^ 1 . In Fig j2] one can observe the fairly good agreement between the full result and this first 
approximation. 
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Figure 2: Numerical evaluation of \G > (x, t) \ at zero temperature T = resulting from bosonization, 
as a function of x and t (solid blue lines). The weight of the sharp peak at x = vpt decreases for 
increasing propagation times t. The dashed blue line shows an approximation of the Green's function 
(see Eq. (43 1) which yields good qualitative agreement with the full solution. The plot is done for 
U q = f7 e- (9/9c)2 with ^ = 2ira = 5 
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xq c 



Figure 3: a) The electron Green's function \G > (x, e)|, for different lengths i as a function of energy 
e, evaluated using bosonization, for U(q) — Uoe~( q / qc ^ [solid blue lines] and for U(q) = [Toe - ' 9 / 90 ' 
[dashed red lines]. At high energies we show the result coming from the semiclassical approach 
for the Gaussian potential [solid black line]. The red area indicates the regime of validity for 
the semiclassical (high-energy) approximation. The interaction strengths are: 2na = 2 (a) and 
2na = 15 (b). In the high-energy limit, the semiclassical solution is valid for arbitrary coupling 
strength. However, the energies for which the description is valid become larger for increasing 
coupling strength. In (b) this limit is beyond the presented energy interval. Therefore here we do 
not show the semiclassical solution. 



3.4 Green's function vs. position and energy 



As shown above in Eq. (15) and (20 1, the current through the interferometer is determined by the 
propagators G > ^ < (e, x). Therefore, in the following our main focus will be on this function, which 
can be thought of as the amplitude for an electron of energy e to propagate unperturbed over a 
distance x. 

The function is shown in Fig. [3j where we plot the numerical evaluation of the exact result 
obtained using the bosonization technique. This is done for two values of the coupling strength a 
and for different interaction potentials. There are some main features which can be observed by 
having a brief look at Fig. [3] while for a detailed discussion we refer to the following section. 
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(i) At x = 0, where \G > (x = 0, e)| = z/(e) equals the tunneling density of states, there is a finite 
dip at low energies. This is a static interaction effect. For repulsive interactions it represents the 
suppression of the tunneling density by a factor vf/v, due to the interaction-induced increase of 
the velocity v. At high energies (e 3> vpq c ), the non-interacting density of states is recovered. 

(ii) At any fixed energy e, the Green's function decays with increasing propagation length x. 
The shape of the decay (as a function of x) becomes independent of energy for high energies. In 
contrast, the decay is suppressed for energies below e ~ vfQc, and there is no decay in the limit 
e — > 0. The decay of the GF is equivalent to dephasing (since in our model there are no interbranch 
interactions and correspondingly no vertex corrections). As a consequence, the absence of decay at 
zero energy will lead to perfect visibility at T — 0, V — > 0. 

(iii) At larger x, there are oscillations in the Green's function. These result from the double- 
peak structure in the time-domain, with peaks at x = vt and x — vt. These lead to a beating term 
exp[iwx(w^ 1 — w -1 )] in |G > (x,o;)|. Therefore, the period of oscillations in the energy domain is 
determined by the difference between the bare and the renormalized velocity (see Fig (3]), viz.: 

2ir ., 11 

Se ss — ; with 77= -. (44) 

xrj v v 

3.5 Large coupling constants 

In this section we want to discuss briefly the shape of the Green's function in terms of the coupling 
strength. We emphasize that, once the shape of the interaction potential is given, the only dimen- 
sionlcss parameter left is the coupling constant a — Uo/(2ttvf)- All the other parameters may be 
absorbed into a rescaling of the result. 

In Fig. [4] we show |G > (e,a;)| as a function of energy for various coupling strengths (different 
curves), both at x — and at some finite propagation distance x ^ 0. For small coupling a > 0, 
we just observe the suppression of the tunneling density of states discussed above. Upon increasing 
the coupling strength, a series of rounded steps emerges, suppressing the tunneling density even 
further. The same features can be seen in the shape of the GF at finite x, though there they are 
superimposed by the decay (describing decoherence) and the oscillations as a function of energy 
(discussed in the preceding section). To identifiy the oscillations in energy which we observe even for 
small coupling strength in Fig|4j; we divide the GF for x ^ showed in Fig|4]D by the tunnel density 
of states. As expected, those oscillations are robust against a change in the coupling strength. We 
have not found any simple analytical model to discuss the structures observed here. However, note 
that in Fig|3j3 one observes that the step structure is more pronounced for the Gaussian potential 
compared to the results for the exponential shape. That shows the strong influence of the shape of 
the interaction potential on the step structure. 

We note that the previous discussion in the (e, fc)-space (as opposed to (e, x)) had found non- 
analytic structures for the case of a box-shape potential U q [M] . 

4 Semiclassical model of dephasing 

Up to now we have applied the bosonization technique in order to get exact information about the 
decoherence the electron suffers while passing through the interferometer. However we actually do 
not know what is going on in more physical terms. First of all, we do not know in detail how to 
distinguish between the single electron we are considering while traveling through the interferometer 
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□ Q Q 

| G > (e, x)\v F | G> (e, x) \v F \G> (e, z) |/i/(e) 




g c ^F 9c^f g c t>F 

Figure 4: Green's function at zero temperature T = as a function of energy e for various coupling 
strengths a, at a; = (a) and xq c = 2 (b). 

c) \G > (e,x)\ for a;g c = 2, divided by the tunnel density of states v{e) = |G > (e,a; = 0)|. This might 
be interpreted as the electron's coherence as a function of propagation energy and distance. 
Here the potential is U q — 2navFe~'-- q ^ q ^ , where the various values of a are (from top to bottom): 
a = 0.16,0.9,1.6,2.3,3 . 
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and the bath electrons which are present at the same time. In a Fermi liquid there is no question 
about the nature of the single electron, i.e. in the vicinity of the Fermi edge it can be described 
as a quasiparticlc. In contrast, in the bosonization approach everything is described in terms of 
collective bosonic excitations, i.e. density fluctuations. Therefore, as soon as the electron tunnels 
into the interacting system, there is no way to trace this special electron any longer. This fact leads 
to difficulties in understanding the decoherence intuitively, since in the moment of tunneling the 
coherent phase information is encoded into the bosonic degrees of freedom. Fortunately, it turns 
out that indeed it is possible to find simple physical pictures which are helpful in understanding 
the process of dephasing in more detail besides the mathematical solution. 

For energies much higher than the Fermi energy (to be made more precise below), it turns 
out that it is possible to describe the interaction of a single propagating electron with all the 
other electrons, by viewing them as a bosonic quantum bath. To see this, in section [4] we apply 
an intuitive, semiclassical method which is able to reproduce exactly the results from the full 
bosonization formalism in the high energy limit. 

In a recent work |36] we have already briefly reported on universal dephasing for high-energy 
electrons at long distances for T = 0, based on the semiclassical approach to be discussed in more 
detail in the following. 



4.1 Semiclassical approach to the Green's function 

Electrons at high energies e propagate at the unperturbed speed vp, as can be observed from the 
corresponding limiting behaviour of the plasmonic dispersion relation. The decoherence, i.e. the 
Green's function, can be described in a transparent semiclassical framework, that becomes exact 
in the limit of high energies. We will confirm this later by comparing against the full bosonization 
solution. Within this semiclassical picture, one thinks of the electron propagating through the 
channel, while accumulating an additional phase due to the interaction with the bath formed by all 
the other electrons. To model the effective, bosonic bath acting on the single electron, we make use 
of the plasmonic dispersion relation which was derived using the full bosonization technique (Eq. 



33 ) . As the electron we consider is flying high above the Fermi sea, we can neglect the backaction 
of the electron onto the bath. In this picture the electron only experiences the intrinsic fluctuations 
of the bath. The potential acting on such a single high-energy electron is obtained by convoluting 
the density fluctuations with the interaction potential: 

V(t)= [ dx'U{x' -vpt)p{x',t). (45) 



Note that this definition implies, that the effective potential fluctuations V(t) experienced by the 
single electron are just the fluctuations of the bath evaluated at the classical electron position 
x = Vpt at time t. This is why we call this model "semiclassical". 

If we were dealing with a classical fluctuating potential V(t), the electron would simply pick 
up a random phase ip(t) = — J dt'V(t'). In that case the non-interacting Green's function would 
have to be multiplied by a factor ^e^ 4 ') to obtain the correct Green's function. However, if the 
quantum nature of the bath becomes important one has take care of the non-commutativity of the 
operator V(t) at different times. This can be done by introducing a time-ordering symbol: 



,-F(t) = 



= ( Texp 



i I dt'V{t') 
'o 
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exp 



dh / dt 2 (TV(h)V(t 2 ) 



(46) 



The time t = x/vp in Eq. (46 1 is determined by the propagation length. This is actually identical to 
the decay of coherence of a single level whose energy fluctuates. In various contexts, this is known 
as the "independent boson model" [3^1 HD] , or the case of "pure dephasing" in a (longitudinally 
coupled) spin-boson model [T]. 

We note that the same kind of approach to dephasing of ballistically propagating electrons has 
been introduced previously, both for a situation with a general quantum bath |171 1211 110] , as well 
as for two coupled Luttinger liquids [30] . 

Furthermore we note that the decay is independent of energy e. This is because the propa- 
gation speed is energy-independent, and the distance to the Fermi edge becomes unimportant at 
high energies as well. Qualitatively, we have seen this feature before in our discussion of the full 
bosonization solution. 

In summary, the decay of coherence, described by F(t), is completely determined by the fluctu- 
ation spectrum \VV) — J dte wt lV{t)V{Q)j of the potential seen by the electron in the moving 



frame. To proceed further we express the time-ordered correlator ( TV{ti)V{t 2 ) ) as a sum of 
commutator and anti-commutator part: 

1 



(fv^vih)) = 2 



{V(h),V(t 2 )} ) + sgnfo - t 2 ) [Vih), V(t 2 )] 



(47) 



The real part of F(t) and therefore the decay of the Green's function depends on the symmetrised 
part of the correlator. This part is formally similar to the correlator of classical noise, though it 
also contains the zero-point fluctuations of the plasmon field: 



Re[f(i)] 



dh / dt 2 ({V(h),V(t 2 )} 



du sin 2 (wa;/2ui?) 

2^ W 2 " 



{V,V} 



(48) 



In addition, a phase — InxF(£) shows up in the exponent. It is due to the commutator of V , and 
thus it represents a purely quantum mechanical contribution. In terms of the Fourier transform of 
the spectrum, this yields 



lm[F(t)} = -Uo dt i Jo d *2 - *a) ( [V(h), V(t 2 )] 



l r+°° dui 

" 2 J-oo 2-7T 



[V,V] 



From Eq. ( 45 ) , we obtain for the potential spectrum in the co-moving frame 



(49) 



(50) 



The argument uj + vfQ indicates that we are dealing with the Galileo-transformed spectrum of the 
density fluctuations. As a result, the spectrum of the density fluctuations gets tilted compared to 
the original dispersion relation (see Fig.[5|. Making use of the plasmonic dispersion u)(q), which we 
obtained from the bosonization method, the density-density correlator yields 



1 - e-P" 



[Q{ q )8{u-u>{\q\))-Q{-q)8(u + u{\q\))] 



(51) 
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We obtain 



^v,v]) u = f§\u^{[p,p]) qtU+qVF 

[{V, V)) u = l£ \U q ? coth ([p, p]) 

The symmetrized correlator can be written as 



q,u)-\-VFq 



{v,v}) = 



d( l ITT 12 ^ ( /M«)\ X f \ U ql\ 

-\U q \ gcoth — )S{ U - — 



Plugging in these correlators in Eq. (48) and Eq. f49[ ), finally yields F(t) 
Re[F(t)\ 



dq coth ( Ml) 



o q 



(U q q. 
1 — cos I — — t 
zir 



lm[F(t)] = f 
Jo 



dq 



2tt 



t - 



dq 2^ 



t . 



(52) 
(53) 

(54) 

(55) 
(56) 



Indeed, up to an additional energy renormalization — J °° dq^, the resulting Green's function 
G > (x = vpt,t) = e~ F ( t 'GQ (x = VFt,t) is identical to the one stemming from the bosonization 
technique for x — vpt [compare to Eqs. (40) and (41 1] . This constant energy shift had been 
incorporated into a redefinition of the chemical potential, as noted in the section on bosonization. 

To be precise, Eqs. ( 55|56 1 show that the decoherence of electrons moving along the trajectory 
x = vpt can be described exactly in terms of the semiclassical approach (which is based on the 
assumption of high-energetic electrons). In order to prove the statement that the decoherence at 
large energies can be understood within the semiclassical description, one has to confirm the initial 
assumption that high energetic electrons actually move with the bare Fermi velocity vp. Already, 
we observed that the Green's function G > (x,t = ^) determines the high-energy 



3.3 



m section 

regime of the Fourier transform G > (x, e). This followed from the analysis of the Green's function 
G > (x 1 t). Beside a broad peak moving with v, one observes a sharp peak in time, moving along the 
trajectory x = vpt, which is obviously responsible for the contributions to G > (e, x) at large energies. 
In addition, we confirm the assumption numerically. As we are not able to perform the Fourier 
transformation analytically, instead Fig. [H] shows the numerical evaluation of the function e~ F ^ 
and the direct Fourier transform of the Green's function G > {x,t) resulting from the bosonization 
in the limit e — > oo. Obviously, those turn out to be identical. Therefore we conclude that the 
semiclassical approach is becoming exact in the limit of high-energy electrons. We emphasize that, 
apart from the transparent interpretation of the decoherence of high energetic electrons via the 
semiclassical method, even the observation that the Fourier transform of the Green's function at 
high energies turns out to be identical to [G y (x,t) /Gq (x,t)]GQ (e, x) with x — vpt is non-trivial 
and serves as a starting point for the analytical discussion of the decoherence in this regime. 

Having established the connection between the semiclassical approach and the full bosonization 
solution, we now turn to the properties of the symmetrized potential spectrum. For this, first 
we focus on the zero temperature case T — (see Fig. [5]). At high frequencies, we obtain a 
singularity ({V, V}}^ oc 1/ -\/w max — \lo\ at the cutoff frequency w max = max(w(g) — vpq), which 
is the maximum frequency in the Galilei-transformed plasmon dispersion relation. Such a maximum 
frequency arises due to the momentum cut-off in the dispersion relation (which results from the 
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Figure 5: a) Plot of the effective spectrum (VV) J' =0 for a > of the plasmonic bath at the particle 
position which moves with velocity vp. The spectrum is linear in u) for small frequencies and 
diverges like ^ 1 _j when approaching w max ■ The inset shows the dispersion relation of the 
plasmonic bath in the laboratory frame of reference as well as in the electron frame of reference, 
where it is getting tilted, [dashed curve] 

b) G > (e,x) as a function of x for large energies e 3> q c VF ■ The solid orange line denotes the 
numerical evaluation of the bosonization result, while the solid black line shows the semiclassical 
result. For a better comparison of the results, the inset shows a blow-up of the oscillations. The 
small deviations are due to finite numerical precision. In addition, the dashed line denotes the 
solution of the Keldysh pertubation theory. All the plots are done for U q — U^e~^ q / qc ^ with 
Si = 27ra = 5. 
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finite range of the interactions) . Due to this smooth momentum cut-off, the velocity of plasmons in 
the limit of large momenta is identical to the Fermi velocity vf of the electron. By transforming the 
potential to the moving frame, it gets tilted (see Fig{5|. As the velocity of the plasmons in the limit of 
large momenta is identical to the velocity of the moving frame, the effective dispersion relation shows 
a maximum. The singularity in the spectrum arises due to the fact that ui(q) m cu niax + uj"(q — q*) 2 /2 
in the vicinity of q* , where co(q*) — w max . Note that an interaction potential U q with a non- 
monotonous decay in q may give rise to several such singularities, corresponding to the local maxima 
of u>(q) — vpq. 

At low frequencies u> <C vpq c , the spectrum increases linearly in uj, corresponding to "Ohmic" 
noise, which is ubiquitous in various other physical contexts pQ. For interaction potentials that 
are smooth in real space (i.e. where all the moments of \U q \ are finite), we find that the leading 



low-frequency behaviour is determined solely by the contribution to Eq. ( 50 1 stemming from small 
q. The result is (here for a > 0): 

({V,V})^ = j^- 2 ^ = 2n\u J \. (57) 
(v — vf) ^ 

Most remarkably, the dimensionless prefactor (the slope) of the noise spectrum turns out to be 
completely independent of the coupling strength a, which drops out. This is in contrast to the 
typically studied non-chiral Luttinger liquids, where an Ohmic spectrum has been found with an 
interaction-dependent prefactor [3D]. As a direct consequence, the electron's Green's function shows 
a universal power-law decay at long distances, as we will discuss in more detail in the next section. 

4.2 Universal dephasing for high-energy electrons 



We insert ( 57 1 into the long-distance limit of Eq. ( |48[ ) : 



rr,/ m n f duj sin 2 ( lux /2v F ) i(~ \' J 
Re Fa; = 2 lim / 1 / F, { \V,V ) 

r c , 1 „ r max dojsm 2 (ujx/2v F ) i(~ T ^1\ T =° 

w const + / duo- +2 — ^ (\V,V\) 58 

Jv F /x u 2ir uj 2 \L )/ u 
v ' 

\u(xlj c /vf) 

where the choice of to c <C w max is arbitrary (but related to the constant), and |G > (e — > oo, x)\vf = 
e -Rc[F(a;)] _ ^ g & con sequence of the logarithmic contribution to the exponent, the leading asymptotic 
behaviour of the Green's function is a power-law decay with an interaction-independent exponent 
1: 

|G > (e,x)| cx -. (59) 
1 1 x 



We recall the fact that Eq. (591 refers to the Green's function in energy-coordinate space, as a 
function of propagation distance x at large energies. Having shown previously that the current and 
thus the visibility only depend on this function (see Eqs. ( 15|20 1), we may denote it the 'coherence' 



of the electron. Therefore, Eq. (591 shows that the coherence itself displays a universal power- 
law decay. Please observe that in the absence of interactions \G > (e 7 x)\ would be constant, so the 
decay indeed implies interaction-induced decoherence. The surprising feature is the independence 
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Figure 6: Coherence of an electron propagating at high energies in an interacting chiral system, as 
a function of propagation distance for various values of interaction strength a: 2ira — 0.1 (solid 
black line), 2ira = 1.0 (solid green line) and 2ira — 5.0 (solid red and blue). The potential is 
taken as U q = Uoe~\ q / qc \ and temperature is zero except for the red line where T/q c vp = 0.005. 
The non-interacting case would give vf \G > (e, x)\ = 1. The long-distance decay is universally 
given by oc 1/x, independent of interaction strength. Note that for decreasing coupling strength 
the asymptotic power-law decay sets in for increasingly larger propagation distances. At finite 
temperatures, this power-law decay turns into an exponential decay for large x with a decay rate 
T v depending on interaction strength (inset). 
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of the exponent from the strength of the interactions (while the prefactor, not displayed here, would 
indeed depend on the details of the potential) . 



In writing down Eq. ( 58 ) , we have neglected the contributions of large momenta in Eq. ( 50 1 . 
These will lead to a subleading correction to the power-law, which we discuss below. The third 
term in Eq. ([58] ) is responsible for oscillations in the coherence \G > (e, x)\, on top of the decay (see 
below). Fig. [6| shows the decay of G > (e — > oo, a;) for different coupling constants a. 

In order to understand how this generic result for the asymptotic decay is compatible with the 
non-interacting limit (a = 0, where |G > (e,a;)| is constant), we have to discuss the range of validity 

of the asymptotic behaviour. As the linear slope in the effective spectrum ^{^i^}^ applies 

only at \u)\ <C w max , we must certainly require u) max x/vF S> 1. Since u max vanishes linearly with a, 
the limiting regime is reached at ever larger values of x when the interaction strength is reduced. 
This shift of the range of validity is shown in Fig. [6] 

The oscillatory modulation is due to the square root singularity at u> — > a-> max in 
Its amplitude depends on the interaction strength |a| but vanishes at long distances (see appendix 



ATI 



cv x , , , / / s n sm(u mim x/vF + vr/4) 

F(x) « const + m(u! c x/VFj — C ■ (60) 

y/2n\a\q c x 

where C is a numerical constant, which solely depends on the explicit form of the potential. 

We now discuss the deviations from the leading low-frequency behaviour of ^{^: ' taking 



into account the contributions resulting from large q in Eq. ( 50 1 , for smooth potentials like U q = 
Uoe^^ q ^ Qc ^ s the correction is ({^'^}) = 27uj/(s ln(|a| VFq c /u)) (A.2I. In real space this 

expression translates into a correction ^ ln(ln(|a|<7 c x)) to the decay function F(t = x/vp) (A. 3 1. 

At finite temperature T ^ 0, the long-time limit is given by an exponential decay \G > (x 1 e)| cx 
exp[— T^x/vf], with a decay rate 



T<p = nT 



1-^ 



= nT\l + a" 1 T 1 . (61) 



This follows from the long-time limit of Eq. (W8k together with Eq. ( 54 1 



RelFix)] ~2r d q ™ a (glW(^» coth (H^0)\) . (62) 
Jo Q V 2T ) 

Using the identity lim^oo Sln = tt5(x), we get Re[-F(:z;)] = ttT\\ + a _1 | _1 a;, from which we 
obtain the aforementioned decay rate r^. In Fig. [6] the decay rate is shown as a function of the 
coupling a. For small a, this rate vanishes as T v — TrT\a\, i.e. it is non-analytic in Uq oc a. Such 
dephasing rates proportional to T have also been found in non-chiral Luttinger liquids |31l 1301 13"2"] . 
At large repulsive coupling, Uq — > +oo, we have the universal result T v — > -kT. For attractive 
interaction, T v diverges upon approaching the instability at a — > — 1, where v — > and where the 
resulting low-frequency modes are thermally strongly excited. 

We note that this behaviour is somewhat surprising when compared to other problems of de- 
phasing. When considering pure dephasing of a two-level system by an Ohmic bath, a power-law 
decay at T — 0, with an exponent 7 set by the coupling, automatically implies an exponential 
decay at a rate T v = tt"/T at finite temperatures. This follows from the fluctuation-dissipation 
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theorem (FDT) which turns the T = Ohmic spectrum into a white-noise spectrum with a weight 
proportional to T. 

However, in the present case, we have to take into account the Galileo transformation, which 
turns the laboratory-frame temperature T into an effective temperature T e g in the frame moving 
along with the particle at speed vf- In order to establish the FDT for the effective potential 
fluctuations in the electron frame of reference we therefore have to define the effective temperature 
T c g by demanding for u> J, 



coth 



V, V 



(63) 



In the low-frequency limit we get ^{^'^}^ = • coth (u) 1 1 + or 1 1 /(2T)) . Together with 

Eq. (52 1, this yields the effective temperature T c g = T|l + a _1 | . In the moving frame the 
frequencies are reduced by q c VF and therefore the effective temperature is also smaller. Only for 
large repulsive interactions (v 3> vp), the transformation does not matter. Therefore, in this limit 
T c g — T and the universal, coupling-independent power-law for T = turns into a universal decay 
rate at finite temperatures. 

Finally, turn briefly to the question of observing these features in experiments. We note that the 
universal power-law decay of the Green's function should be observable in principle in the Mach- 
Zchndcr interferometer setup, as it directly translates into a decay of the visibility itself. To obtain 
a numerical estimate, we assume a screened Coulomb potential with screening length q^ 1 ~ 10~ 7 m 
and a finite channel width b ~ 10~ 7 to. Following [25J the edge state velocity can be assumed to be: 
vf ~ (10 4 — 10 5 ) — . Then the dimensionless coupling constant a is of the order 1. To reach the 
high energy limit, the applied bias voltage V has to fulfill: q e V q c VF- For vp = 10 4 — , it turns 
out that one has to apply V ~ 10 2 [iV which is in the range of the bias voltages typically applied 
in experiments (in [T2] V ~ lO 1 /^)- We mention that the long-distance limit should be reached 
for x WF/^max ~ 10 _1 /iTO, which is shorter than the typical interferometer arm length (e.g, in 
[T2~] X12 ~ 10 1 ^m). Thus there is some hope that the high-energy as well as the long-distance 
limit is accessible in the experiment. However, we note that the magnitude of the Green's function 
G > (e, x) gets suppressed strongly when reaching the long-distance limit, i.e. x ff/o; max . For 
a = 1, the magnitude of |G > (e,x)| yields: |G>(e,x)| - 10^' (see Figje}. With help of Eq. Q this 



translates into a visibility of the order vj ~ 10~ 4 (for this see also section [fj]). Therefore, the direct 
measurement of the power-law dependence of the visibility for large bias voltages seems to be a 
challenging task, unless one finds a way to optimize further the relevant parameters. 



5 Keldysh perturbation theory 
5.1 General remarks 

The semiclassical method we introduced in the preceding section provides us with an intuitive 
picture for the case of electrons flying high above the Fermi sea. To complete the analysis, in 
the present section we employ perturbation theory to discuss the behaviour at short propagation 
distances and weak coupling for all energies. This includes the low energy regime, where the 
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Figure 7: The relevant processes which contribute to the self-energy up to second order in the 
coupling. The wiggly lines indicate the interaction. Note that the plasmonic diagram and the 
vertex correction differ only by an additional exchange process. That produces a minus sign, such 
that the diagrams tend to cancel each other at low momenta. 



influence of the Fermi edge becomes important. We will apply Keldysh (i.e. nonequilibrium) 
perturbation theory up to second order in the interaction strength a. The main outcome of the 
perturbation theory is that the tunnel density of states is affected by renormalization effects while 
the decay of the GF in the close vicinity of the Fermi edge is suppressed. We will find that the 
suppression of decoherence is brought about by a cancellation between two second order diagrams 
near the Fermi edge. 

During this subsection we fix the notation, following the review [41] . Using the Keldysh time 
leads to an additional matrix structure of the GF which reflects the fact that one has to differentiate 
between points in time which lie on the backward or the forward branch of the Keldysh contour. 
Beyond this additional structure, all the known Feynman rules remain exactly the same. After 
performing the rotation in Keldysh space |41j . the representation of the matrix GF G and the 
related matrix self-energy is given by 



G = 



' G R 


G K ' 











G A 


£ = 





Y, A _ 



(64) 



where we introduce the Keldysh GF, G K (x, t, x' , f) = -i([^f(x,t),^(x',f)}). First we will derive 
an expression for the retarded self-energy T, R (e.k), which will be used to calculate G R that can 
be related to the single particle propagator G > in equilibrium. Starting from the matrix Dyson 
equation G(e, k) = Gq(c, k) + Go(e, k) ■ S(e, k) ■ G(e, k), one finds that the retarded Green's function 
only depends on the retarded self-energy: 



G R (e,k) = 



eo(k) 



1 



E*(e, k)' 



(65) 



In the following we calculate the diagrams up to second order for a linearized dispersion relation, 
but for finite temperature and for an arbitrary interaction potential. In the end we compare the 
results of the perturbation theory with the results of the bosonization technique. The relevant 
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processes are shown in Fig (7] There are two second order diagrams, which can be identified as 
the interaction with a plasmonic excitation and a corresponding diagram containing an additional 
exchange process (that can be viewed as a vertex correction diagram). The crucial point is that 
the vertex correction counteracts the plasmonic processes in the vicinity of the Fermi edge, leading 
to a suppression of the decay of the GF. 



5.2 Evaluation of the diagrams 

The starting point of the calculation is the evaluation of the the unperturbed electronic propagator 
matrix Go- In addition to the usual retarded and advanced GFs of free electrons, G (&,k) — 
UJ - V pk±m+ ' tne Keklysh propagator is given by G(f(LJ,k) = — 27ritanh(4^) ■ S(uj — vpk). In con- 
trast to the advanced and retarded GFs, the Keldysh propagator contains information about the 
electronic spectrum as well as about the occupation of those states. Therefore at this point one 
could introduce arbitrary non-equilibrium states which is the main advantage of working on the 
Keldysh contour. However, as we are describing channels which are only weakly tunnel-coupled to 
each other, we will calculate equilibrium Green's functions. 

As in second order we explicitly include interactions with free plasmons, we also derive their 
propagators here. To this end we identify the bosonic field with the potential V(x, t) = J dx' U(x — 
x')p(x', t). This is identical to the potential we introduced in the semiclassical description [although 
in the latter case we specialized to x = vpt]. The bosonic propagators are defined as D R / A (x,t) — 
TiO(±t)([V(x,t),V(0,Q)]) and D K (x,t) = -i({V(x,t),V(0,0)}), respectively. A straightforward 
calculation yields for the retarded and advanced propagators 

y ,q > q J 2TT(LU + e Q (k)-e (k + q))±iO+ q 2tt y ,q > y ' 

As we assume the system to be in equilibrium, we use the FDT to obtain the plasmonic Keldysh 
propagator: 

D K (uj 1 q) = 2icoth(f3u/2)Im[D R (u,q)}. (67) 

Now we can proceed calculating the various contributions to the self-energy up to 0{a 2 ). Con- 
sidering all the possible Feynman diagrams, we are left with the first order Hartree-Fock diagrams 
and in second order with the plasmon diagram and the vertex correction (see Figj7]). Thus we can 
express the self-energy as = Snartrcc + Sp ck + Spi asmon + Syertox- Those can be evaluated 

according to the rules given in [JT]. We use the equal-time interaction propagators U R / A (q) = U q 
and set U K = 0, as usual. 



First-order contributions: Hartree-Fock diagrams The Hartree diagram yields a global 
energy renormalization Ai?Hartroo = U(q = 0)p. In the framework of the Luttinger model the 
electron density diverges as there is no lower boundary of the electron spectrum. However, formally 



one can include the energy shift into the definition of the chemical potential (see subsection 3.1l. 
In the following we omit the Hartree contribution. 

The self-energy contribution due to the Fock diagrams is 

E«( w> fc)=-i- J ' dqU q f(k-q), (68) 
which for zero temperature yields ^p T = (k) = —^U(x = 0) + ^ dqU q [see Figj^i]. 
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The resulting fc-dependent energy shift describes, in particular, the renormalization of the eletron 
velocity near the Fermi edge. This also affects the tunneling density of states, leading to a suppres- 
sion (for repulsive interactions, a > 0) or enhancement (a < 0). This can be seen in Fig{3j Again, 
the constant shift can be incorporated in the definition of the chemical potential. 



Second-order contributions: Plasmonic excitations and vertex correction The electron's 
coherence decays by interacting with the plasmons, i.e. the density fluctuations of the other elec- 
trons. The plasmon diagram (see Fig. [7]) represents one of the contributions describing this physics. 
It yields 



E?(e,fc) 



(dq) / (dw) [G«(e - w, k - q) ■ D«{u>, q) + G*(e -u,k-q)- D*(u, q)] , (69) 



where the second term contains the Fermi function, which introduces the effects of the Fermi edge 
on the coherence. Inserting the propagators, the contribution can be written in a compact form 



8tt 2 



dq 



[coth(/3/w F <j/2) + tanh(/?/w F (fc - q)/2)} 



(70) 



which for T = reduces to S 



P (e, k) = G«(e, k) ■ ^ /l fe| dqVjq. Thus, at T 
vanishes for k — > 0. We note in passing that the structure "coth + tanh" generically occurs in 
discussions of dephasing, where it describes both the strength of the thermal fluctuations and the 
influence of the Fermi function, i.e. the physics of Pauli blocking |42l [T3 [2TJ [lOj HI |9] . In the Umit 
of high energies, the result of Eq. (70 1 can be rewritten in terms of the potential fluctuations at 
the particle position, as discussed in the preceding section. Specifically, we have lim^oo Ep T=0 = 

Gq ■ (V(x = 0,t= 0) 2 ). For a plot of the function ££ T=0 , see Figj8t. 

Finally we derive the vertex correction, mentioned above, which alter a rather lengthy calculation 
yields (see appendix [B| 



this contribution 



Eg 



(e,fc) 



G§(e,k) 



16tt 2 



dqi 



dq 2 U qi U q , 



tanh( ^ (fc ~ gl ~ 92) )) 



tanh( 



/3v F ((k-qi), 



-tanh( ^ ((fc ' g2) )) 



■tanh( ^ (( ;~ g2) ) tanh( ^ F(( ;~ 9l) )-l 



(71) 



This expression simplifies for T — to E^ r=0 = — • ^2 Jq^ dqi /iju_ gi dl2U qi U q2 . The con- 
tribution from the vertex correction as well as the total second order correction to the self-energy 
(jp_i_Y = Ep Ey are shown in Figj8] The crucial feature is that up to second order in momentum 
k the plasmon diagram and the vertex correction cancel exactly against each other, whereas for 
high momenta [k 3> q c ] only the plasmon contribution remains while the vertex correction tends 
to zero. In summary, when calculating up to second order in the coupling, we already see that 
the dephasing is suppressed in the vicinity of the Fermi edge. As the comparison with the exact 
bosonization solution shows, this conclusion holds true qualitatively to all orders. We can the pre- 
ceding expressions to evaluate the retarded Green's function G R and from this the propagator G > . 
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Figure 8: (a) The GF from bosonization (solid blue lines) vs. the results from second-order Keldysh 
pertubation theory (dashed red lines) for different propagation distances as a function of energy. 
Note the fairly good agreement, (b) Fock contribution to the self-energy S^(fc), which starts with 
a linear slope (corresponding to the electron velocity renormalization) and then saturates for large 
momenta, (c) Total second order contribution the self-energy (dashed, black line); plasmonic (solid 
black line) and vertex correction (lower, red line) separately. The cancellation for small momenta 
is discussed in detail in the text (in all figures U q — U a exp(—(q/q c ) 2 ) and 2na = U /v F = 2) 

The final results has evaluated numerically. The results are illustrated in FigjHJ Expanding the 
Green's function for small propagation distance x <C vf (^(O) 2 )^ , yields: 

\G>(e - oo, x)\v P « 1 - i<y(0) 2 )( — ) 2 + ■■■, (72) 

2 vp 

which coincides with the expanded exact result. The good agreement of the bosonization result 
and the Keldysh perturbation theory for small |or| is shown in Fig. [8] as well as in Fig. [5] However, 
for large x the perturbation theory fails. 

For a detailed study of chiral interacting electrons of the spectrum in (e — fc)-space, starting 
from the bosonization result, we refer the reader to |34j . 

5.3 Summary of the Keldysh perturbation theory 

To summarize, the perturbation theory shows that there are two different energy regimes: In 
general, the GF decays as a function of propagation distance due to the interaction with the density 
fluctuations. This is particularly pronounced at high energies e 3> vpq ci where we have also shown 
that the Keldysh result and the semiclassical (or bosonization) approach coincide at short distances. 
For low energies e <C VFq c , the decay is suppressed. On the other hand, at low energies the Fermi 
velocity is renormalized due to virtual processes, leading to a modification of the tunneling density 
of states. It is important to note that these two different energy regimes only emerge since we are 
dealing with interaction potentials of finite range. 
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Figure 9: a) The amplitude of the coherent (flux-dependent) part of the current / co h through the 
interferometer for different arm lengths x. The curve for x = is equal to the flux-independent 
current / C oh(^ = 0, V) = Io(V) (red line), which implies that the interference contrast is perfect at 
zero armlength. 

b) The visibility vj as a function of bias voltage V for various arm-lengths x\ = xi — x. The red 
line denotes the semiclassical calculation. The small deviations from the bosonization result vanish 
completly for larger V. The plot is done for U q — U$e ^^i/i^' with ^ = 2ira = 3. 

6 Visibility and Current 

As mentioned in the beginning, the results for the GF G > (e,x) we worked out in the foregoing 
section can be applied directly to the evaluation of the current and the visibility (Fig{9j). 

Fig|9^i shows the current through the interferometer as a function of voltage for different arm 
lengths. Here we restrict the considerations to the symmetric case, i.e. Xi = Xi- For x — the 
coherent part of the current obviously is identical the flux-independent part, which implies a perfect 
visibility (at T = V = 0). The suppression of the current at small voltages is due to the velocity 
renormalization which lowers the tunnel density (for repulsive interactions). However, as the change 
in the tunnel density influences the classical and the coherent part in the same way, it does not 
show up in the visibility at all, i.e. vi(V, x = Q) = 1 (see Figj9|j). 

The behaviour of the coherence of the GF for high energies discussed above (using the semi- 
classical approach) can be transferred directly to the discussion of the visibility. Therefore, in the 
limit of high voltages V and T = the visibility is determined by the factor |G > (e — » oo, x)\ 2 . This 
follows from the fact that for higher voltages the contribution of the high-energy electrons becomes 
dominant. It also implies that the visibility at high voltages becomes voltage-independent. Corre- 
spondingly, the exponential decay at finite temperature is transferred to the visibility vi as well. 
For i^O the dephasing reduces the coherent (flux-dependent) part of the current which leads to a 
decrease of the visibility. At small voltages |<7 e V| <C q c ^F the visibility decays only very slowly with 
increasing interferometer length (see the discussion of the Green's function). In the limit V — > 
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the visibility is approaching unity vj — * 1 , which is consistent with the fact that in equilibrium and 
at zero temperature there is no dephasing. In contrast to the Green's function itself, which shows 
oscillations as a function of e, the visibility does not show pronounced oscillations as a function of 
V. 

7 Conclusions 

In the present paper, we have studied dephasing by electron-electron interactions in a ballistic 
interferometer. We have considered the case of electrons moving inside ballistic, one-dimensional, 
chiral channels, such as edge channels in the integer quantum Hall effect. The interference contrast, 
in the limit of low transmission at the beam splitters, can be expressed via the Green's function 
of the interacting system. We have studied the decay of the electron's coherence as a function 
of propagation distance, employing three different approaches: The exact bosonization solution, 
a semiclassical approach that becomes exact for high energies, and diagrammatic nonequilibrium 
(Kcldysh) perturbation theory. Our most important physical result is that at high energies the 
decay of coherence at T = becomes a power-law with a universal exponent, independent of 
interaction strength. We have also shown that second-order perturbation theory compares well 
with the exact solution at modest coupling strength and short distances. This may be important 
for potential applications to more elaborate setups that cannot be treated exactly any more. 
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A Semiclassical approach 

A.l Amplitude of the oscillations in G > (e,x) as a function of x for T = 



The third term in Eq.(58l is responsible for the oscillations in G > {e — > oo,x), while the main 
contribution stems from the square-root singularity at uj = \uj max \. This contribution yields: 



{*•*}). 



s|uwl 



const. — Re 



(dw) 



e "p 



|«w|-a 



(73) 



From Eq. ( 54 1 we obtain the w-dependence of the correlator in the vicinity of the singularity 

U 2 (q )qo 



where to(q*) = ma,x(\u> q — v F q\), q = q* + y 



iax| - W) 

t = Id 2 /l^lg^ 

5 — 2 dg 2 ^ 27T i 



result into Eq. ( 73 1 yields 



9=9 



(74) 

Plugging in this 



1 U\q*)q* 



S7T 2 w 2 

max 



Re 



|^ max | 



e 



1 C/ 2 (<f)g* 



max 



• Re 



|ow|-<5 



(75) 



Evaluating the integral in the limit x — > oo gives 

1 tf 2 (<f)<? 



5F(x) 



S7T 2 w 2 



Sinf W m(U t h 



V F 



tt/4) 



sin(|w max |^ +tt/4) 

"O ■ 



(76) 



where C denotes a numerical prcfactor which depends only on the form of the interaction potential, 
not on the interaction strength a. 



A. 2 Sub-leading corrections for smooth potentials U q at T = 



The subleading corrections result from the contributions of large q in Eq. ( 54 1 (here u> > 0) . We 
start with the expression 



{v,v} 



dq U 2 qS(u} 



q\u q \ 

2ir 



)• 



(77) 
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The sub-leading correction to the low-frequency behaviour of {{V, V}) for a smooth potential 
V q = U e^^ q / q ^ s is given by 



{v,v}) 



(sub) 



dq U:q 



S (l - Isub) 



f q (\U q \q) 



(78) 



where g su b fulfills: 

<7sub 



_ q^b £ \q sab /q c \ _ Therefore, as q 3> q c for uj small enough we obtain for 



<?sub = Qc In 
Finally, the subleading corrections yields 

<{*■*}>; 

for w <C |a|wF9c- 



2ttlu 



1/3 



(sub) 



2ttuj 



(79) 



(80) 



A. 3 Sub-leading correction in time-domain for T = 

The subleading correction to the decay function F(x) results from the expression = |^ — ^ x J( 2vf ^ . 

f * f. ~| \ (sub) J I 

< V, V > ) . With Eq. (80) (omitting an additional constant) it yields 



(du) 



2-kuj 



(81) 



where for reasons of brevity we set |a|uF9c = w*. As we can choose the cutoff frequency oj c such 
that u c /ui* < 1, we can assume < u> < 1. This enables the substitution co = e~ x : 

F {suh \t) = / dx- = 

s J \ti(u>*x/vf) x 



^ rln(u}"x/v F ) ^ ^ 

- / dx — — - ln(ln(|a|g c a;)) — const. 

s JlniW/uc) x s 



(82) 



B Keldysh perturbation theory: Vertex correction 



The vertex correction to the retarded self-energy has the structure (compare to the review of [41): 
'A) 2 j^yi J dh J dk 2 J d^ j duj 2 U qi U q , x r 



Ey(e,fc) - \ ~: —j i ,} J "-(J 11^92 

G^(e - W2, fc - fc 2 )G^(e - wi - w 2 , A - *i - fc 2 )G^(e - w x , A - fci) 

+G^e - w 2 , fc - fc 2 )G^(e - - u; 2 , fc - fci - k 2 )G%(e - uji, k - fci) 
+G^'(e -uj 2 ,k- fc 2 )G^(e - w x - w 2 , fc - fci - k 2 )G$(e -w^k-fa) 
+Gg(e -u 2 ,k- k 2 )Gfi(e - w x - w 2 , k - fc x - k 2 )G^(e — Ui,k— fci)] 



(83) 
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